Exact solutions for the periodic Anderson model in 2D: A 
Non-Fermi liquid state in normal phase. 

Peter Gurin and Zsolt Gulacsi 

Department of Theoretical Physics, University of Debrecen, H-4010 Debrecen, Hungary 

(January, 2001) 



Abstract 

Presenting exact solutions for the two dimensional periodic Anderson model 

with finite and nonzero on-site interaction U > 0, we are describing a rigorous 

non- Fermi liquid phase in normal phase and 2D. This new state emerges in 

multi-band interacting Fermi systems above half filling, being generated by 

a flat band effect. The momentum distribution function together with its 

derivatives of any order is continuous. The state possesses a well defined Fermi 

energy (ep), but the Fermi momentum concept is not definable, so the Fermi 

surface in /c-space is missing. The state emerges in the vicinity of a Mott 

insulating phase when lattice distortions are present, is highly degenerated 

and paramagnetic. A gap is present at high U in the density of low lying 

excitations. During low lying excitations, quasi-particles are not created above 

the Fermi level, only the number of particles at ep increases. 
PACS No. 05.30.Fk, 67.40.Db, 71.10.-w, 71.10.Hf, 71.10.Pm 
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I. INTRODUCTION 



The topic of non-Fermi liquid (NFL) behavior in D > 1 dimensions and normal (non- 
symmetry broken) phase (NP) is currently of great interest!. This is mainly due to the 
large amount of experimental results, obtained in principle in the last decade, showing NFL 
behavior in the NP of a variety of materials, including D > 1 dimensional systems of 
great interest. Examples are: high T c superconductors!, heavy- fermionsi, layered systems! , 
quasi-one dimensional conductors, doped semiconductors, systems with impurities, materials 
presenting proximity to metal-insulator transitions!!, etc. These results changed considerably 
our understanding of interacting Fermi systems. Indeed, until recently, Fermi-liquid (FL) 
theory seemed universally applicable to all sufficiently pure interacting Fermi systems, and its 
main features even to dirty systems, provided that their NP is not destroyed by a symmetry 
breaking process^. This "dogma" has been based on high precision experimental verifications 
in liquid He3 and simple metalsll. The concept of FL itself has been introduced by Landau 
many decades ag (for a thorough discussion seel), and in principle has the meaning that 
in spite of the interactions, the low energy behavior can be well described within a picture 
of almost noninteracting quasi-particles. Formulated in rigorous termsSJi, in a normal FL 
we have a one-to-one correspondence between the non-interacting and interacting single- 
particle states (determined e.g., by a perturbation theory convergent up to infinite orders). 
Furthermore, a quasi-particle pole is present in the single-particle propagator that gives rise 
to a step-like discontinuity of the momentum distribution function at the Fermi surface, 
whose position is specified by a sharp Fermi momentum value kp. The observation of NFL 
behavior in the materials presented above polarized a huge intellectual effort in the last 
decadeS for the understanding of this new fermionic state. In this field the theoretical 
interpretations are often based on multi-band mo delsE, the presence of a some kind of gap 
in the NP being clearly established in many cases and subject of intensive experiment alEl 
and theoreticaflEl studies. However, despite the great number of papers published in the 
field (see for example the references cited in! or!), and the fact that the observed most 
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interesting and important normal NFL properties emerge in two spatial dimensions (2D), 
(for example the normal phase of the high T c superconductors), on the theoretical side, for 
pure systems, the existence of a NFL state in a NP has been proved exactly only in one 
dimensions (i.e. Luttinger liquid^). The extension possibility of NFL-NP properties to 2D 
has not been demonstrated rigorously up today. In fact, a rigorous theory of a NFL normal 
state in higher than one dimensions is missing. 

Driven by these state of facts, we started to focus our attention on possible NFL states 
using exact methods which are applicable in higher than one dimensions as well@. Based on 
the obtained results and a conjecture made by us@, in this paper we are reporting the first 
rigorously derived NFL state in 2D. We deduced for this reason exact solutions for a real 
space version of the periodic Anderson model (PAM) in 2D. The model is analyzed on a two 
dimensional square lattice, in case of non- vanishing and finite Hubbard on-site repulsion U . 
In the paper in fact two qualitatively different solutions are described: a completely localized 
and a non-localized one, which represents the first exact results reported for PAM in 2D and 
finite U. The solutions are valid on two surfaces of the parameter space, i.e. on restricted, 
but continuous and infinite regions of the T = phase diagram, extended from the low U 
to the high U regions up to U = oo at U > 0. 

The derived non-Fermi liquid state is given by a flat band effect in multi-band systems 
with more than half filling. The obtained properties are extremely peculiar: the system in 
case of the described solution possesses a well defined Fermi energy ep in conditions in which 
the kp Fermi momentum cannot be defined, and the momentum distribution function is 
continuous together with its derivatives of any order. The system has also a natural built 
in gap, which is clearly present in the density of low lying excitations at high U. The state 
is paramagnetic and non-insulating. The gap symmetry is a possible symmetry allowed by 
the described 2D lattice, and depends on the starting parameters of the system. The state 
emerges in the proximity of a Mott insulating phase when lattice distortions in the unit cell 
are present. During low lying excitations quasi-particles are not created above the Fermi 
level, only the number of particles increases at e^. 
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Concerning the flat band features (FBF), we mention that such characteristics have been 
clearly observed in different systems where strong electron interactions and strong correlation 
effects play a main role. On numerical side, FBF are present for example in results connected 
to 2D Hubbard mode@~El, or 2D t — J modefi Experimentally FBF are seen in angle- 
resolved photo-emission (ARPES) data of high T c cuprateSH For layered systems ARPES 
often shows main bands without any sharp characteristics in n^, or give results interpreted 
via FBF assumptions^. Band structure calculations for these systems often reflects a Fermi 
level positioned exactly at the bottom of a conduction band with large effective mass around 
its minimum, below which a gap is present!!. We further wish to mention that connections 
between superconductivity and FBF were also clearly pointed out by Imada et al.@, and FBF 
can be seen as well in experiments related to heavy-fermion materials^. On the technological 
side, for example Lammert et alMlhave shown that squashing carbon nanotubes, FBF can be 
achieved around ep, where a mismatch of nearly isoenergetic k states may have unexpected 
application possibilities. 

The remaining part of the paper is constructed as follows: Section II. presents in detail the 
analyzed model and the general form of the deduced ground-state wave-functions. Section 
III. characterizes the obtained solutions from the point of view of the k - space representation 
of the Hamiltonian and wave- vectors described, and Section IV. analyzes magnetic properties 
of the system in the studied ground-states. Section V. describes a completely localized 
insulating solution, Section VI. presents the new non-Fermi liquid state in normal phase 
and 2D, Section VII. summarizes the paper, and the Appendix, containing mathematical 
details, closes the presentation. 

II. THE MODEL AND GROUND STATES DEDUCED 

We are describing in this Section the model we use and the ground-states detected for it 
in restricted domains of the phase diagram. 
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A. The model 



We are considering in this paper a 2D square lattice described by a two-band model whose 
Hamiltonian for the start is given in direct space, containing on-site repulsive interaction in 
one band. The starting point will be sufficiently general in order to give us the possibility 
to characterize in detail the state we are presenting. However, the model contains also 
restrictions. Based on physical considerations, and denoting by d n the positions of the nth 
neighbors of a given but arbitrary lattice site, we are taking into consideration in this paper 
only the n < 2 (i.e. nearest, and next-nearest neighbor) contributions in the Hamiltonian. 
With these considerations, our starting Hamiltonian can be given as 

H = H + U , H = f c + f f + V + V + E f , (1) 

where the non-interacting terms have been denoted together by H . With i denoting an 
arbitrary 2D lattice site position fi, the contributing terms in Eq.fll]) can be explicitly written 
as follows 

Tc = E (^c]^ £ j+x,<x + h.c.) + E ( *c, y 4<r a J+y^ + h -c.) + 

E ( tc,x+y Cj,<r Cj+(x+y), CT + h.C.) + E ( t c ,y-x Cj+x.fr Aj+y.a + h.C.) , (2) 



J>f 3," 



Tf = E (WjV W + h.c.) + E + ^ c -) + 

E ( */,*+!/ /j+(y+x),<7 + /i.e.) + E ( */.»-* /jW /j+y,<x + /i.e. ) , (3) 



J,°- J,o- 



ft = £ ( d l W + ) + E ( /2a a j+x,. + ) + 

E ( Kl 4. W + ) + E ( 5. *W + h - c - ) > ( 4 ) 



^2 = E ( C L /j+(x+y),<x + h.c. ) + E ( V lx+y fL £ j+(x+y),<x + h.c ) + 

E ( ^2,y-x C j+x,(j /j+y,o- + ^-C ) + E ( ^2,y-x fj+x,a C j+y,cr + h.C ) , (5) 



= E 04,4. + h.c.) , E f = E f J2flh«, U = Uj2M^n(l- ( 6 ) 

j,cr j,cr i 

In these expressions x (y) represent elementary displacements along a single bond of length 
equal with the lattice constant a in x (y) direction, respectively. Taking into account different 
couplings along different lattice directions, we allow in fact the study of the system with 
distorted unit cell as well. The interaction term U is the usual Hubbard interaction, where 
h[ a = /; a is the particle number operator for f i a electrons. 

Our starting H from Eq.(|l|) represents a prototype of an interacting two band system 
in 2D. It describes a square lattice containing fermions distributed in two bands b = c, f. 
The kinetic energy term is in fact f = J2n=i J2{ n ;ij),b,a i.^bd n Kah,° + ^- c -) > w here t b j are 
hopping amplitudes, a is the spin index, and (n; i , j ) has the meaning of a sum over bonds 
connecting nth neighbors, every bond being taken into account once. The hybridization 
between the bands is composed from on-site Vq and nth neighboring sites hybridizations 
V — En=i,2 Kt- As can be seen, in Eq.(P the non-local hybridization V = V\ + V2 
contains the nearest-neighbor (Vi) and the next-nearest-neighbor (V^) contributions only, 
the V n , n > 3 long-range terms being considered negligibly small. The on-site energy Ef 
for the second band fixes the relative position of the two bands, and finally, U represents in 
fact the on-site Coulomb repulsion in the second band b = f, making the model interacting 
(U > ). We mention that U is present on all lattice sites. In order to help the reader in a 
better understanding of the notations, we are presenting in Fig.l. the hopping amplitudes for 
c-electrons (A), and hybridization matrix elements for c\fj type transfers (B). The hopping 
amplitudes for /-electrons, and the hybridization matrix elements for f}tj type transfers are 
similar. For them, only the index c has to be changed in / in Fig.l. A., and the superscript cf 
has to be replaced with fc in Fig.l.B. As given above, our Hamiltonian H can be considered 
as a 2D - PAM given in a direct space version, or a 2D two-band Hubbard model containing 
the contribution of the Hubbard U only in one band, the other band being non-interacting. 

We concentrate in this paper on a specific region of the phase diagram of the model 
Vjj. This region Vfj can be defined as the domain of the parameter space where H can 
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be expressed through elementary plaquette contributions in the form A\ a Aj a . Here, the 
plaquette operator Aj a is build up from a linear combination of the starting b annihilation 
fermionic operators acting on the elementary plaquette /. We are denoting the four positions 
available on a plaquette by j = 1, 2, 3, 4 , starting from the down-left corner and counting 
anti-clockwise (see the notation of sites from Fig.l.). The Jth plaquette operator for a fixed 
spin a can be generally expressed as Aj )(7 = J2j,b a j(/),6^j(/),o-> where, because of the translation 
invariance, djini — a j,b plaquette independence will be considered for the emerging eight 
coefficients present in A/ jCr , that have to be deduced. Taking the translational invariance 
explicitly into account, the plaquette operator acting on four corners of an elementary square 
plaquette / becomes 



Ai, a = axj fi + a 2J f 2 + a 3)f f 3 + a 4J f 4 
+ ai,c ci + a 2 , c c 2 + a 3iC c 3 + a 4iC c 4 . 



(7) 



Starting from Eq.(|7|), and taking into consideration periodic boundary conditions, the prod- 
uct A\ a Aj a summed over plaquettes for a fixed spin index a can be written as 



E A-I,a A I,o 

I 

Na 



E 

i=l 



N A 



H,c\ 



E4. 



E { [ ( a tc a 2,c + al,c «3, c ) c\ p c j+Xj(T + h.c. 



+ E k/ 

V i=i 

+ 



NA * 

E/ fj,cr fj,v 

j=l 



[ ( a *,c «4, c + a 2 c a 3:C ) c\ a c j+yi(7 + h.c. 



+ 



( <A,f a 2,f + a\j a 3J ) / j+X)CT + h.c. 
( a ij a 4j + a*2j «3,/ ) lie /j+y.a + h - c - 

a l,c a 3,c cj,(T ^j+(x+y),(J + h.C. 
a lJ a 3,/ /j jCr /j+(x+y), CT + h.C. 

( a l,c a 2j + al,c a 3 ,f ) c\ a / j+X)CT + h.c. 
( a* j a 2fi + a\j a 3 c ) fj a Cj +Xj(J + h.c. 

( < c a 4,/ + a2, c «3J ) 4<x /j+y^ + hx - 



+ 



+ 
+ 

G 2,c a 4,c Cj jCT Cj_|_(y_ x ) j0 - + /l.C. 
a 2J a 4,/ /j+(y-x), CT + h.C. 

+ 
+ 
+ 



+ 



+ 
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( a[ j a^ c + a* 2 j a 3>c J fj a Cj +yi<J + h.c. 



+ 



<c «3J 4<r /j+(x+y),<r + h.C. 
a *2,c «4J c[ a / j+ ( y -x),<x + /l-C. 



+ 



l l,/ «3 iC /j )(T Cj + ( x+y ) iCT + h.c. 



a 2,f a 4:,c /j l<7 Cj+(y-x),cr + fr-C 



+ 



^ £ a tc 04J ) f £ 4" ^ j ' CT j + ^- C " | 



(8) 



In Eq.(||), A^a denotes the number of lattice sites. A comparison of Eq.(|T]) and Eq.(§) shows 
that H can be expressed via the A I)C operators as H = — YLi a A\ a A It(7 + K N, if the 
following conditions are satisfied 



tc,x 




0-4 )C «3,c , 


- r C)2/ — a l c a^ c -+- a 2 c a3 >c , 






a *.f a 2,f + 


a\j a 3 j , 


_ = a lJ fl 4,/ + Ct2,/ 03,/ , 




^c,x+y 


= a* l c a3 jC , 


tc,y—x ~ 


= «2,c a 4,c , — tf,x+y — a lJ a 3,f i 


- tf,y-x - a 2J a 4J , 


- v,% = 


- a lc d2j + 


a l,c a 3J ) 


— Vifj. = a* j d2 )C + d\f a 3,c ) 




- < - 


-- a* lc a 4 j + 


a 2,c fl 3,/ j 


— V/ y = a\ j a^ c + j a^ tC , 




V 2,x+y 


= a l,c a 3J 


- 1/ /C 






V 2,y-x 


— a 2,c °4,/ 


V 2,y-x 


— a 2 j o-4 )C , 




-V = 


4 

2J a tc a iJ ' 
i=l 


4 

^ = E 

i=l 


1 I 2 


(9) 



E 



f 



(10) 



i=l 



In the expression of -ffo the operator A" = J2i,b,a represents the total particle number 
operator. Based on Eq.(^), it can be observed that A\ a A I a . + A I (7 A\ a = Yh=i [ \ a i,c\ 2 + 
- 1 = 2 AT — Ef , so the Hamiltonian ffn can be written as 



H,f\ 



H = £ 44 + ifiV - 2[2iT - ^] AT A 



(11) 



In the presence of the interaction, besides H Q from Eq. (|TT]) , we have to take into con- 
sideration the Hubbard term U as well. However, we may observe that U can be exactly 
transformed as 



S 



u = up' + u C - ^^a, (12) 

i,cr 

where P' = P( an d P/ = ( 1 — ^i,2,T — ^1,2,4. + Wi,2,T^i,2,l) ■ m ^ e decomposition presented 
in Eq.fll2f) the P' operator is a positive semidefinite operator. The reason for this is simple. 
Pj applied to a wave function gives one if on the i site there are no / electrons present, and 
gives zero, if on the site i there is at least one / electron present. As a consequence, P' 
representing a sum of non-negative numbers, it is a positive semidefinite operator. 

We have further to observe that in Eq.(^) the U M,a term simply renormalizes the 
Ef contribution from the Hamiltonian. Keeping this information, introducing the notations 
Ef = Ef + U , and G = J2i,a A/> ? the starting Hamiltonian from Eq.fll]) can be 
written as 

H = [G + UP'] - [U N A + 2(2 K - E f ) N A - K N]. (13) 

The decomposition presented in Eq.(|T3"D is valid if the conditions presented in Eq.@ are 
satisfied, and in Eq.(^) the Ef value is replaced by Ef, where 

Ef = K - Y, \Hf?- (14) 

i=l 

We mention that from mathematical point of view, the parameter space region is given 
by the solutions of the system of equations Eqs.(|9|,P^). If this system of equations admits 
solutions for the coefficients a^b from Eq.(^), we are situated inside Vfy. 

From Eq.(|l3D can be seen that in conditions presented by Eqs.([9lJl4j), the analyzed Hamil- 
tonian from Eq.(|l]) can be written as 

H = P + Eq , (15) 

where the positive semidefinite operator P and the constant number Eq are given by P = 
G + U P' , and 

E% = KN + ( U + 2 E f - 4K)N A , (16) 

where N being the eigenvalue of N, represents the number of particles within the system. 
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B. The detected ground-states 



For the Hamiltonian from Eq. flT5|) , P being a positive semidefinite operator, the ground- 
state energy is = Eq , and the ground-state wave function is that | *$?g ), for which we 
have P | tyg ) = 0. The interesting aspect from the mathematical side of the problem is 
that | ) can be explicitly expressed in the concentration range situated at and above 
3/4 filling, obtained in this way a many-body, fully quantum mechanical solution in the 
interacting case and 2D. For this, we have to take into consideration that starting from the 
definition presented in Eq.(|7p, for plaquette operators we have A\^ a A\, a , + A\, a , A\ = . 
Furthermore, the number of elementary plaquettes / from the system is equal with the 
number of lattice sites, and the product = Yii m / based on the operator (see also 

Eq.©) 

hf = («i,T At + «u hi ) ( 17 ) 

containing arbitrary constants, creates an / electron on every site of the lattice. As a 
consequences, for example at 3/4 filling (i.e. iV = 3iV A ), the ground-state (not normalized) 
wave function | ) becomes the ordered product 

l*?> = A F^\0), (18) 

(3=1 

where 

f« = n 4t . f {2) = n A ii . (is) 
/ i 

and |0) being the bare vacuum with no fermions present. Indeed, because of the A\ a A\ a — 
property, we have G\^ v g ) = given by the F^F^ product in Eq. (|T%D , while F^ introducing 
on every lattice site one electron, obliges | ^ ) to have at least one / particle on each site, 
preserving p> | \&u ) = o as explained below Eq.(H). 

We would like to underline here that the a^ a coefficients in Eq.(^7|) are completely 
arbitrary, so the most general form of the ground-state is in fact a linear combination of 
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components present in Eq.([H|) taken over all possible combinations of all possible a- h(7 values. 
This means that in the most general case, the F- h f operator (see Eq.(jT7|)) contained in 
has the form 

hf = E 9{a-^} ( «i,T At + «u hi ) > (2°) 

{«i, CT } 

where g{ aia } represents numerical coefficients (connected to fixed {cti i(7 } sets) restricted only 
by the normalization to unity of the whole wave function. 

The solution for the ground-state can be also written for the system doped above 3/4 
filling with arbitrary 1 < n r < N\ number of electrons. For this reason we must define 
the operator 

^i (4) = e T ( i? T + e** 4, T ) + £l(fZ + ^ 4 d ) , (21) 
which allows us to introduce randomly n r additional electrons in the sistem above 3 N\ via 

= E C {1} fl ^ (4) • (22) 

{i} I 

In Eq. (p2|) , represents an ordered product containing n r arbitrary chosen lattice sites 
taken as a possible combination of n r sites from ]Va possibilities, and C^j are numerical 
coefficients. As a consequence, in the doped case with 1 < n r < A^a, the ground-state wave 
function becomes 

\< d ) = A ^|0>- (23) 

(5=1 

We mention that along this paper we will restrict our study to the N > 3 N\ (i.e. n r > 0) 
case. 

It is extremely important to mention that the presented ground-state wave functions in 
Eqs. ( |18| , |23| ) are valid only in interacting (U > 0) case, and cannot be perturbatively 
obtained from the U — noninteracting limit. The reason for this is that at U = 0, as 
explained below Eq.tJiTf), the P|\l/^ =0 ) = property is entirely given by the F^F^ product, 
the F^ operator being completely arbitrary. In the noninteracting case the operator P 
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reduces to G, so we obtain inside Vfy at U — the equality G \ ) = with the 
ground-state wave function 



/3=1 



Q\0), (24) 



where Q is an arbitrary operator. As it can be seen, the concrete expression of the F (3) 
operator is determined exactly by the nonzero U > value of the interaction (see below 
Eq.(|19])). Based on these characteristics mentioned, we would like to underline that is no 
way to re-obtain (together with all expectation values that it gives) | tyg ) from | ^/° g ) (or 
vice versa) in the U — > limit. 

The remaining part of the paper is devoted to the study of the physical properties of the 
| ) wave function. We stress that depending on different solutions allowed by the system 
of equations Eqs.([||l4|), the ground-state wave function given mathematically in Eq.(|i~8l) 



describes from physical point of view even qualitatively different ground-states, which will 
be analyzed below. 

III. THE if-SPACE REPRESENTATIONS 

First of all, we would like to understand the physical background of the Hamiltonian form 
presented in Eq. (Tf3|) . We remember, that H from Eq.([T3|) represents an exact representation 
of the starting Hamiltonian from Eq.(^) in conditions in which Eqs.([| |T^) are valid (and 
admit solutions). 

A. The Fourier transform for H 

In order to answer this question, let us transform the starting Hamiltonian from Eq.(|]) 
in fc-space. In order to do this systematically, let us first concentrate on Hq. Denoting by 
k the two-dimensional reciprocal space vector, and using for the operators the b = c, / 
notation together with the Fourier sum b- 3 = £V exp — ik • rj , the kinetic energy 
terms becomes Th = Y.? e b r bt b? , where 
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(V^e"'*^ + % +2/ e +i ^+^) + (t^e-^-M + t* bty _ x e +a <^ ) . (25) 

Concerning the notations, we mention that using d n introduced at the beginning of Sec. II., 
defining t b3n = \ t b ^ \ exp(i (/) tbli ), we simply have e£ ff = 2 I t b j n I cos^bjL - £• d n ). 

For the on-site energy at the / level we simply obtain Ef = Ef J2k a f% a fka- Introducing 
the notation T' f — TV + Ef together with the definitions e% — e% , e£ — Ef + el , 

J ■> k,a k,<T k,a ■> k,a 

we simply have at the level of Fourier transforms T c + T'f = J2 b k a e \ a t ( no ^ e ^ ne 
difference between el and el ). 

k,cr k,a' 

In the case of the hybridization, the on-site term becomes V Q = Eg c ( ct f% a + h.c. ) . 
For the non-local hybridization V, as shown in Sec. II., we have V = V\ + V2, where, for 
n = 1, 2 we obtain y„ = Ej£ CT ^fco- + ^- c - ) • The hybridization matrix elements 

in these expressions for n = 1,2, are given by 



v l,k = 


( ^ e" 


- i k-x 


V 2 ,k = 




e -ik 


( V 2,y-x 


e -ik-(y- 


x) _|_ 



".V 



J 



'2,j/-a: 



/c * e + . (26) 



Introducing the notation = Vo + V l % + V 2 t, the total hybridization can be given as 
V + V = ( V % cl /, + h.c. ) , and for H we get 



^0 = E [ C + ^ f u + ( V S 4 ct + fc. c . ) ] . (27) 

k,a 

The diagonalization of Hq can now be simply done. For this reason we have to introduce 
the row vector W$ = ( ct , /J ) , and the (2 x 2) matrix M with components Mi 1 = 

rC,(X ft" ^ CT rC,CT 

e|, Mi j2 = V^, M 2i 2 = ep M 2j i = V^. The H from Eq. ([27|) will contain under the sum over 
k the expression ( W\ M Wk ) ■ The diagonalization of H in fc-space reduces to the secular 
equation written for the matrix M. We obtain from this ( e| — E^ (el — J — | V^\ 2 = 0, 
from where, as expected, two bands arise (i = 1, 2) 

1 



fc,i 2 



+ (28) 
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The expression of in Eq. fl28|) is given by Tjg = y ( e| — el ) 2 + 4 | \ 2 . If we are using 
now Eqs. (|9],[10]) in the expressions given for e|, el, and in terms of the starting Hamiltonian 
parameters in Eqs.(^,^), we realize that the following equality holds 

\V S \ 2 = {K - el){K - e[). (29) 

This is a lengthy but straightforward calculation which in fact, easily can be done. Using 
now Eq.([29"l) for T^, we find T| = ( e| + el — 2 K) 2 , so the band structure given by Eq.([28|) 
becomes 

E %1 = K = constant , E l2 = e| + el - K . (30) 

From Eq.(||) we can see that K > 0, since K is positive definite. Furthermore, we have 
K > e| as well. Indeed, using Eqs. (p|,p5|) we obtain 

4 

X - e| = 53 I a i,c I 2 + [(a*, c a 2, c + al c a 3iC ) e~ * kx + c.c. ] + [ ( a* )C a 4 , c + a 3 ,c ) x 
i=i 

e _<fc » + c.c] + [a^ c a 3iC e- i( ^ + ^ ) + c.c. ] + [ a^ c a 4 , c e~ i ( " k * > + c.c.]. (31) 
Introducing now the notations 



. ■ k X + ky . . fc^ — fc^ 

/ + J 



a l,c = Q l,c e 2 > a 2,c = a 2,c e 2 , 

_ . k x + fcy . fc x — k y 

a 3,c = a 3,c e 2 } a 4 c = a 4iC e + 2 ; (32) 
the exponential factors disappear from Eq.([JT]), and for K — e| we find (see also Eq. ([!(])) 

iv - = ( a lc + a 2 c + a 3>c + a 4iC ) ( a 1>c + a 2jC + a 3)C + a 4>c ) = 

I a 'l,c + a 2,c + a 3,c + «4, C I 2 • ( 33 ) 

However, K > e| via Eq. (|29|) means > as well, since | V% | 2 is a non-negative 
number. As a consequence, the band structure obtained in Eq.(30) contains an upper 



band that is completely flat ( 1 = E\ = K ) , and a lower, normal, k - dependent 
band ( E% 2 ) with dispersion. The k dependent gap between these two bands is given by 
Ajf = Ei- E^ 2 = 2K - (e c n + el ) , and based on Eqs.(||, 0) > holds. Introducing 
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A = Min[ At], we have A > 0. Being important below, at this point we mention that 

Ag > means A > as well. From physical point of view, A > means that the 
diagonalized bands from Eq. ([3(]) are never intersecting. In case of A = , Et t and Et 2 
touch each other, while for A > the diagonalized bands are completely separated. 

From Eq.(|5DD it can be observed, that the conditions from Eqs. (P|jrO|) that allows the 
transformation of the starting Hamiltonian Eq. ([!]), into H from Eq.(|l|), are exactly the 
conditions that give a band structure containing a completely flat (i.e. dispersion-less) 
upper band seen in Eq.(|3C|). 

When the system is interacting and the U > Hubbard term is present, the Hamiltonian 
becomes H given in Eq. fll3|) . In the ground-state, because of P' \ ^ ) = 0, the effective 
Hamiltonian has in fact exactly the form of H^, excepting that Ef is renormalized as 
Ef = Ef + U (i.e. the condition from Eq.([T0|) has to be changed to that given in Eq.(Jf4|)), and 
the energy scale, as seen from Eq.([R]), is shifted with U N\. As a consequence, effectuating 
the band structure calculation as presented above, using instead of Ef the Ef value, we 
re-obtain (shifted with U N\) for the ground-state the structure presented in Eq. fl3"0|) . The 
same holds for excited states which give P' \ ^ ) = as well. 

It is interesting to mention at this point, that a two - band system with a band structure 
as given in Eq.(^), above half filling (i.e. with more than two electrons per lattice site), 
has a well defined Fermi energy positioned at ep — Ei — K (where K is a k independent 
constant), but the Fermi momentum kp is not definable. A such type of system has no 
Fermi surface in fc-space, so for this case, the Luttinger theorem is without meaning. 

Another aspect that has to be accentuately underlined, is the fact that the Hubbard 
interaction gives effectively its contribution in the flattening of the upper diagonalized band 
Ej: x in the case in which the here reported solutions are valid. In order to understand this, 
first let as mention that for U = 0, as seen from Eq.(^), the solutions from Eqs.( |T8| , p3D are 
not applicable. Let turn then to the U ^ case, and analyze a concrete pedagogical example. 
Consider for example the Hamiltonian parameters t CjX = t CtV = 12, t C)X+y = t CtV _ x = 
— 4, tf^x = tf, y = 3, tf yX+y = tf^ y - x = —1, Vo = 18, V\ >x = V\ w = —6, V2, x +y = ^2, y -x = 2, and 
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Ef = 2. The written numbers are expressed in |i/ )X +jf| units. This particular case describe 
an isotropic situation, with V c f = V^ c = V, which has been chosen only to be easy 
for the reader to follow, without diminishing the general physical content of the behavior 
reflected by Eq. fl30|) . The Et^ band from Eq.(|28|) with U = and parameters given above, 
is presented in the upper plot of Fig. 2. (the k x and k y values cover the first Brillouin-zone, 
i.e. [— 7r, +7r] in units of the lattice constant). As can be seen, for this non-interacting case 
we have a normal k dependent band, the system is metallic, and the system of equations 
Eqs.(|9|,^0|) has no solutions (i.e. at U = we are situated outside of Vfy). Hovewer, turning 
the interaction on (note that for U ^ we have Ef instead of Ef in Eq . fl28|) , and Eq . (|i~0"|) 
has to be changed with Eq.(|14|)), the presence of the interaction starts to flatten the E^ 1 
band. Indeed, as seen from Fig. 2. (middle plot: U = 10, bottom: U = 20), the E^ 1 
surface starts to flatten with increasing U values. The completely flat Et l case given in 
Eq.(pOl) is obtained for U = 25, which at the level of the system of equations Eqs.(j9||14|), is 
represented by the solution a liC = 3 + y/E, a 2iC = —2, a 3>c = 3 — y/E, a^ c = —2, a\j = 
-(V5 + 3)/2, a 2,f — 1, ^3,/ — (V5-3)/2, a 4 j = 1. When this solution emerges at 
U = 25, the Hubbard interaction has pushed the system inside Vfr, so its role is fully active 
in obtaining the interacting ground-states described here. 

B. Decomposition into composite operators 

We have to mention that obtaining the diagonalized band picture for H presented in 
Eq.(|3~0~D, the mathematical description can be given in fc-space in term of new (composite) 
fermionic operators, which creates composite fermions in the upper and lower band. To see 
this, we note that the following relation holds 

e l4 c % +elft h + (iVUe^ct fs + h.c.) = 

k k,cr k,<7 k J k,cr J k,<7 \l K l k,cr J k,(7 J 

KCt 6r, + (ef + el - K)Cl Cr, , (34) 

k,l,a k,X,er \ k,a k,a ' k,2,a k,2,a ' V / 

where for j = 1,2 the operators CV . are defined as 
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d U« = 7^ [jK^k* ~ (-lYe+W* ^4%,) ^ • (35) 

From Eq.(pl|) it can be seen that 4>v represents the argumentum of Vk, i.e. Vk = | Vk | e l ^ v . 
Furthermore, 9\ and 82 represent two arbitrary, independent, and k independent phases. 
Using Eqs. (|3"0| , |3~3] ), the operator Ho can be written now in a simple form 

^0 = E ( % 3e,i„ + E t2 C\ v ) . (36) 

The expression presented in Eq.(|3~6"D shows that the operators defined in Eq.(|35|) annihi- 
late (their adjoint create) particles from (into) the ,,diagonalized bands" j = 1 (upper flat 
band), and j = 2 (lower band with dispersion), respectively. Since j C^^ a , Ct , j = 

%' hr <W' > { ' } = { ' ^Lv } = ' the 3fej> °P erators represent 

new, rigorous, canonical and anti-commuting Fermionic operators. 

The transformation into Eq. (|36|) physically is more deep as seems to be at first view. In 

order to see this, we have to analyze the Fourier transform of the plaquette operators Ai )(T . 

Starting from Eq.(^), the Fourier transform of A\ a can be written as 

= E e* fc> [ X c (k) c\ a + X,{%) fl ] , (37) 

k 

where, for b = c, / we have 

X b (k) = al b + a\ h e iU + a* 6 e^ 5+ ^ + a\ h e^ . (38) 

Effectuating the product over the ordered plaquette index I in Hi A\ a (such a product 
emerges in the ground-state wave function from Eqs.(|l8|,[^|)), given by the anti-commutation 
rules of the starting /, c operators, only a product containing different k indices survives. 
Introducing the notation Y = J2p ( — l) p exp[i (r*i ■ kn + r 2 ■ k i2 + ... + rW A • km A ) ] , where 
J2p denotes a sum over all possible permutations of (1, 2, N\) to (zi, i 2 , in a ), and p 
represents the number of pair permutations in a given P, we obtain 

II 4. = Y II [X c $) c\ a + X f (k) flj. (39) 

I k 
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As can be seen from Eg. (|39"D , the contribution into the norm of the ground state wave function 
given by JJj A\ a is proportional to 1Z = H? IZz where IZz = [ \ Xjk) | 2 + | Xf(k) | 2 ] . 
However, using Eq.(^) together with the definitions given in Sec. III. for e| and Eqs.(p|,|38|), 
we find 



X c (k) 



K 



X f (k)\ 2 — K — el, X c (k)X}(k) 



(40) 



For example, in the case of the first relation from Eq.(^), we simply must use X*(k) from 
Eq.(|38l) and multiply it with exp[(i/2) ( k x + k y ) ] in order to see the expression of (K — e|) 
presented in Eq.(p3|). The remaining two equalities can be similarly deduced. 

From the first two relations of Eq. ([40]) we observe that the contribution of the operator 
r y ' or PW) into the norm of the ground state wave function is determined via Eq.(p9[) by 



K % = A S = \X c (k)\ 2 + IX;^)! 2 = 2K - (e c n + el). 



(41) 



Let us now introduce instead of the initial plaquette operators A Ia , new plaquette op- 
erators ,, normalized to unity", i.e. whose contribution under the Yli product give unity into 
the norm of the wave function. In this case instead of Aj a we must use 



B 



Y 



i-i \ x A 



l.a 



Ar a . 



In this case we have instead of Y\j A\ ct , the product 



n si = n 



k L 



X c (k)c\+X f (k) ft 



where 



n Bt , 



(42) 



(43) 



B 



fe,2,<T 



X c {k)c\ a + X f (k)f i 



T 

k,a 



(44) 



Now let us denote X c (k) = \X c (k)\ e^ c , X f (k) = \X f (k)\ e**f , and Vg = |Vg| e i(t > v . Us- 
ing these notations, from Eq.(^) and first two equations from Eq.(^) we find Ct, 



A^ 1/2 ( | X c (k) | c 



- I (p v 



Xf(k) | ft ) . From Eq.(f4"4"|) we have however B 



k,2,a 

t 

k,2.a 



1/2 ( | X c (k) | e* *> ct a + | X f (k) | e* *f fl a ) = e* * A g 1/2 ( | X c (k) | c\ a + \ X f (k) | e* ( */ " *■ ) ) . 
From the third relation of Eq.(|(]) via Eqs.( p9j^0D we obtain e i ^ c ~^^ = — e l<t>v , which 
means e l ( ~ <M = 



e * ^ v , and the expression of Bl becomes 



B 



t 

fc,2,i 



6 °fe,2„ 



(45) 



This means that the Fourier transform of Bj a (i.e. the normalized Az>)> is in fact Cr 2(7 ^ 
we fix the phase 6*2 from Eq.flSo]) to 62 = — <p c - 

e~ l( t> c c\ , which creates 

fc,l,<7 



Similar to Bl , we can introduce the normalized Bl 

k,2,a k,l,a 



a particle in the upper band, by fixing the 9\ phase as 9\ = <p c . We obtain 

^{ X ;Ck)il-X' c (k)fl). (46) 



Since in comparison with C^-, only a phase factor difference emerges, the canonical anti- 
commutation relations remain true for Br . a (J = 1, 2) operators as well. Using Eqs.( [45|Ji"6| ) 
the initial operators can be also expressed via 

_ Mk) + X* f (k) 5j 



(47) 



With the composite operators introduced, for the H Q operator we obtain 

Ho = Y Y E? . Bl Br . , 

and the ground-state wave functions presented in Sec. II. becomes 



TT ( Bl Bl 

11 V fc,2,t fc,2,| 



F® I ) , 

£(3) £(4) J 0) ; 



(48) 



(49) 



where the ordered product over k has to be taken over the whole first Brillouin zone. The 
norm of the ground-state wave functions from Eq.(f49"D is entirely determined by with 
P = 3,4. 
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We however stress, that the k representations for the deduced ground-state wave func- 
tions presented in Eqs.(^) are valid only if the initial ground-state wave functions from 
Eqs.([T^J2^) have nonzero norm. This is the case if rig ^% 7^ holds (see Eq.(f4*2"|)), i.e. 
the diagonalized bands E^. =l2 are completely separated. In order to find such situations, 
mathematically is sufficient to have for all k values from the first Brillouin zone |X c (fc)| > 0. 



IV. MAGNETIC PROPERTIES 

Continuing the study of physical properties of the deduced ground-state, the following 
natural step is to analyze its magnetic properties. This is motivated by the existence of 
ferro-magnetism in some flat-band mo Because of this reason, we have to check if 

something similar is present for our solutions, or not. Starting from this motivation, we 
obtain however paramagnetic properties (i.e. large spin degeneracy) for the deduced ground 
states. The reason for this is the following one. 

As can be seen from Eq . (|2~0"D , the obtained ground-state wave functions have a large 
degeneracy. With this observation in mind, let us consider first the N = 3 N\ (i.e 3/4 
filling) case. We have a d g — 2 Na fold degeneracy in | ^ ) due to the up or down 
orientation possibilities of the iV/3 = electrons from the system, since 2 N\ particles 
fill up completely the lower band, so their contribution in the total spin is zero. From the 
point of view of the a- ly(T coefficients entering in F^> as given in Eq. (|17D , this means that 
it is possible for us to obtain d g linearly independent contributions in | ^ ), by choosing 
arbitrary a^j and a- h i coefficients. The S z value of these states (S being the total spin) is 
situated between [ — N\/2 , + N^/2 ]. The contributions from these states with strictly 
S z = N\/2 can be obtained by choosing = 1 and = for all i sites. This maximal 
S z for the system can be achieved in an unique way, therefore there is an unique state (apart 
from the trivial ( 2 S + 1 ) fold degeneracy) among the linearly independent contributions 
in the ground-state, with S = N\/2. 

After finding this contribution in | ^ ) characterized by S = N^/2 total spin, let us 
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chose one single lattice site ii with a^j = and a- n ^ = 1, the remaining sites i 7^ ii being 
maintained with = 1, a- h i = 0. Moving the position of the ii site along the lattice, 
we obtain all contributions in the ground-state with S z = N\/2 — 1. Let us denote these 
states by | ^p), where (3 — 1, 2, ...N\. These states are linearly independent (see Appendix), 
and they span an dimensional subspace since a single particle spin can be flipped down 
independently on N\ lattice sites. The mentioned | ^ p ) states are also components of 
I tyg ). Since between the | typ ) wave vectors it must be one term corresponding to the 
( S z = A^a/2 — 1, S — N\/2 ) state, we have — 1 linearly independent states with 
S = N\/2 — 1. As a consequence, | ^g) contains not only components with S = but 
also components with S = N^/2 — 1. Continuing this procedure, one can see that every spin 
value S is present among the ground state contributions A^!(2S , +l)/[(A^/2-|-S'-|-l)!(A^/2 — 5")!] 
times (and certainly everyone has 2S + 1 components with different S z values). The S = 
subspace has the greatest degeneracy, i.e. it has the greatest statistical weight. Therefore, 
in the thermodynamic limit the expectation value of the spin goes to zero, i.e. the system 
is paramagnetic. 

In fact, since states with different 5* can be obtained from | \E^) simply by modifying the 
values in Eq. flTT|) , it means that calculating the ground-state expectation value ( S 2 ) 
based on Eqs.(|l^,^), we obtain an expression that depends on the arbitrary coefficients 
a^ a . This means that for | ), the value of the total spin is not fixed, it possesses a 
large spin degeneracy, i.e. is paramagnetic. This property will not be changed even if we 
take into consideration the operator in the ground-state wave function (i.e. doped 
system), because we add in this case in the expression of (S 2 ), contributions depending 
on ot- ha multiplied by arbitrary constants present in Eq.(p2|). As a consequence, the 
detected solutions describes paramagnetic phases. 

Before closing this section, we have to mention that large spin degeneracy characteristics 
of the ground-states for strongly correlated systems have been also reported elsewhere. An 
interesting result on this line, is that described by Arita and AokaS. These authors have 
found for t — t' type Hubbard models ground-states holding simultaneously total spin S = 
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and S = S max values in the thermodynamical limit. In their case, the singlet S = 
component is created by a spiral spin state with a spin correlation length as large as the 
system size, which accompanies the fully polarized ferromagnetic state with total spin S max . 
We have to underline, that our case differs semnificatively from that presented in RefS In 
contrast to Ref.S, in the ground states from Eqs.(IIRBl), all total spin S values are present. 



As will be exemplified later on (see for example a concrete solution described in Eq.fl57)), 
given by the arbitrary nature of the coefficients «i j(7 , a local spin periodicity presence in the 
ground-state is rather accidental, instead of a characteristic (or general) property. 



V. THE INSULATING PHASE 

Besides the fact that the states we obtained are non-magnetic, their physical properties 
still remain to be clarified. These properties depend in fact from the possible solutions of 
Eqs.(f||l4]). To have an insight in these possibilities, let us analyze Eq.(|9|). In this study we 
have to consider the coupling constants from the Hamiltonian H known variables, and to 
try to deduce based on them the parameters entering into Aj t(7 presented in Eq.(^). We 
may start solving this problem by introducing the parameters 



Pi = 


ycf 
*/,aH 


H/ 

-y 


* 

a l,c 

a* ' 
a l,/ 


P2 = 


yfc 
tf,x+y 


d3,c 
a 3J 


Qi = 


ycf 


-X 


a 2,c 


Q2 = 


yfc 
V 2,y-x 


G4,c 


tf,v- 


X 


a 2J 


tf,y-x 


a 4J 



(50) 



Comparing Eqs.([7|,|50D, we can rewrite the plaquette operator A I>a as follows 
4r,<7 = aij ( /i )Cr + p{ c x ,a ) + a-2j ( h, a + q{ c 2 , CT ) + 

03,/ ( H,a + P2 C 3 , CT ) + a 4 J ( /4,cr + ?2 C 4)CT ) . (51) 

The physical nature of the ground-state is strongly influenced by the particle distribution 
created within the lattice by | ^ ) or | ). Since introduces one electron on every 
lattice site, this particle distribution is essentially determined by the product F^ 2 \ But 
on a given lattice site we independently can introduce two electrons with opposite spin, 
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reason for which we focus on the behavior of F^ 1 '. In order to have an image about these 
aspects, we are interested to analyze the situation in which the operator = Yli A] a 
introduces two electrons (one c and one /), at least on a lattice site. If a such situation 
emerges, creates a non-uniform particle distribution within the lattice. This is because 
it introduces in the system N\ particles on N\ lattice sites creating somewhere a double 
occupancy, which must be followed by an empty site. Using Eq.(^l|) and multiplying two 
arbitrary neighboring plaquettes, for a site taking place in both plaquettes the following six 
type of nonzero products may emerge 

h = (fi, a + Pi cL) ' (fl + 11 4<r)> /2 = (fl + Pi ' (fl + P*2 A*) ' 

h = (fl + Pi C) ■ Ul + q* 2 C) > h = Ul + qi c{ a ) ■ (fl + p* cU , 

h = (fl + qi c{a) ■ (fl + q*2 h = ill + Pi ' (fl + <& O ■ (52) 

If at least one Ij term from Eq . (|52|) is nonzero, we identified at least a site i, where II/ A\ a 
introduces two particles. Effectuating the products in Eq. ([^) and taking into consideration 
the anti-commutation rules for the fermionic creation operators, for Ij from Eq.(|52"D we find 

h = Ac fi,a(Pl ~ qi ), h = c\ c fl {pi - p* ), I Z = c{a fi,c ( Pi ~ ?2 ) ) 

h = &l fl ( qi ~ Pi ), h = C fl ( ft - q* 2 ), h = &l fl (Pl~ qi)- (53) 

Usually for the hybridization matrix elements we have = Vj c = V Y for an arbitrary 
element Xi reason for which we focus in the remaining part of the paper to this case (the 
results presented up to the end of Sec. IV. remaining valid for the general Vf* ^ V^ c situation 
as well). This means that we have in fact 

P = Pi = P2 , q = qi = <?2 • (54) 

Now if we consider only real coupling constants in the Hamiltonian H, in Eq.(^) only the 
terms containing 5 = (p — q) survives. From physical point of view, if 5 = 0, the wave 
function | ) introduces on every lattice site the same number of electrons. Since hopping 
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matrix elements for a such type of state are all zero, the ground-state is completely localized, 
i.e. represents a Mott-insulator. 

After this step, considering p and q pure real variables (i.e. all coupling constants from 
H are pure real), we have to check if the solutions allowed by Eqs. (|5|,P^) gives 5 = or not. 
For this reason two cases have to be enumerated. 



A. The symmetric case 

In the first step we analyze the p = q situation, called symmetric case. Physically a such 
type of parameter region can be obtained, for example, when the lattice has a non-distorted 
unit cell. In these circumstances h tX+y = h,y- x and V2, x + y = V%y-x, i-e. (p — q) = 0, 
so all Ij = from Eq. (|53~D . (We underline that rigorously the condition p = q means 
^2,x+ y Vz,y-x — t f,x+yt c , y -x, which exceeds in fact the non-distorted unit cell situation.). We 
mention here, that for p = q the system of equations from Eq. (|9[) admits only solutions of 
the type p = p*. In this case = Ym=i a *ij (p%c + fl,a)i anc ^ the product Y\i from 
the ground-state wave function can be replaced by F a = Hi Fi,a, where = (pc\ a + fl a ). 
We note, that the product in F a is a product over lattice sites instead of plaquettes. The 
contribution into the norm of F a is well defined and equal with (1 +p 2 ) Njv . The ground-state 
wave function becomes for the insulating case 

l*i > = II [hlFiAFij] |0>. (55) 

i 

We mention that this ground-state is present only at 3/4 filling, p = q, and real coupling 
constants in the Hamiltonian. It is interesting to note that for p = q and real, the k 
representations for the ground-state wave functions from Eq.(^) have no meaning since 
(see Eq. ( [4*2| , |63"D ) we have TT^ ^£ = for this case. This means, that the norm of the initial 



| ) wave function from Eq.(|18D is not well defined (i.e. is zero) as expressed in the form 
of a product over plaquette operators in its F^ part. However, a re-structuration in 
the ground-state wave function by changing Aj a to F- XG eliminates from the wave vector the 
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source of the zero norm, and the ground-state from Eq.(|55D becomes well defined, and its 
norm is finite. 

From Eq.fl5"5|) we obtain | ) = Yii | ), where 

F{ = [p 2 c[i 4t fh + P 2 a ^ d ll 4t ftl + P a i.T fll 4t /i!t + P a U d li fi,i fi,i ) ■ ( 56 ) 

Introducing the notation |«i| 2 = (c^f | 2 + la^jj 2 , and q\ = [ ( 1 + p 2 ) |«i| 2 all nonzero 
ground-state one-particle expectation values can be given as 

(flafi,*) = <7ib 2 |ov| 2 + |«i| 2 ], (c\,aCi,a) = [ ^ |«i| 2 + ki^H, 

(4,a fi,<r ) = QiP |tti,-a| 2 • (57) 

Summing up over a we obtain 

E<iU*> = ^> E<^> = ^r> E(4X) = -^r- (58) 

Evidently, here (...) = (*j )/(*/ |*j ) holds. Based on Eq. (||) we re-obtain N = 3iV A , 

and the ground-state expectation values of different Hamiltonian terms becomes 

2pV N A 



(T c ) = 0, (T f ) = 0, (V) = 0, (Vo> , 2 

1 + p 2 J ' 1 + p 2 



1+p 2 

(E f ) = E f N A [ 1 + — 1- ) , (C/> = ff^. (59) 



1 \ [/ A^v 



Summing up all contributions in Eq. (|59|) we obtain for the ground-state energy of the insu- 
lating phase 

H = YT7 2[U + Ef{p2 + 2) + 2pKo] ' (60) 

Introducing (Ri oc ) = ( U ) + ( Ef ) + ( Vo ) as the contribution of the on-site (i.e. localized) 
Hamiltonian terms into the ground state energy, and ( R mov ) = (T c ) + (Tf) + ( V ) 
as the contribution in the ground-state energy of the Hamiltonian terms connected to the 
movement of particles within the system, we find 

( Rioc ) = Eq, ( R mov ) = , (61) 
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which clearly shows that the system is completely localized. 

In fact, from Eq. (fp we have E$/N A = -U + 2 (U + E f ) + [2pV -p 2 (U + E f )]/(1 + p 2 ) , 
which becomes E^/N A = -U + 2 (U + E f ) - K via 

4 

K = p 2 S, V = -pS, U + E f = K - S , S = ^2\ a i,f\ 2 - ( 62 ) 

i=i 

But, for p = q and all coupling constants real, Eqs.fl62|) can be directly obtained from 
the system of equations Eqs.(|,|l4] ) . As can be seen, the ground-state energy obtained in 
Eq.(|60l) gives exactly back the value deduced in Eq.(|i~6l) at iV = 3 N\. Taking into account 
from Eq . (|50|) , that p = V2, x +y/tf,x+y, based on Eq.(|62D, we find the surface (U + Ef)/Vo = 
(1 — p 2 )/p in the {U,Ef, V } parameter space where the described solution emerges. The 
multiple conditions seen in Eq.(|50D relating the hopping matrix elements and neighboring 
hybridizations to the p parameter can be achieved in the simplest way in the case of the 
non-distorted unit cell. As can be seen, 5 = for p = q and p = p* represents a possible 
physical solution for the analyzed problem. 

Before continuing, some observations has to be made here. First of all, from Eq. one 
can see that the sum over a of the nonzero ground-state expectation values presented in 
Eq.(|57p are independent on the a-^ a coefficients introduced by in Eq. (|TTD . This implies 
that the results given in Eq.(^) are valid also for the most general form of the ground-state 
wave function expressed with F^ constructed via Eq.(^) instead of Eq . (|T7|) . We underline 
that the independence on the a- l>a coefficients of the ground-state expectation values summed 
up over a is a general property of the ground-state wave function, and is valid not only in 
the insulating case described here (see Appendix). 

Secondly, we mention that doping the system, the properties described here will be 
destroyed. The reason for this is that introducing electrons above 3/4 filling, the number 
of particles per site will not be constant along the whole lattice, and as a consequence, the 
ground-state expectation values of the kinetic energy terms become to be nonzero. 
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B. The non-symmetric case 



Both p and q being considered pure real, we concentrate now on the possibility of a p q 
solution. We arrive to this inequality for example taking into account distorted unit cell, 
which gives usually the p ^ q condition. The 6 value introduced after Eq.(f>Il) is nonzero, 
so the state we are analyzing is clearly a non-localized state. Solving however for this case 
Eq.(|[), and calculating the norm of the ground state wave function | ) from Eq . fll8|) , we 
obtain zero value. 

Indeed, solving the system of equations Eqs.(|9|JT0D for p ^ q and both real (see Eq.(|33"D), 
we finally obtain 

I *c\ 2 = K - 4 = 4p 2 \t Lx+y \ ( C^j + ^ C^s ) 2 , 
I X f \ 2 = K - e[ = A \t f)X+y \ ( C W + r a f , 

S Q T 

v k = ~ A P ( Ci, s ,fc + — C 2jSj fc ) • ( C l s k + r s C 2)8 j ) , 

K = 2p 2 \t Lx+y \ + V = -2p\t Lx+y \ ^ 1 + r 2 ^ J , (63) 

where s = tf )X /tf jV , r = tf }X /(2tf >x+y ), and, for s > we have C is ^ = cos[f •(£+(— l) t+1 y)] 
and for s < 0, C is ^ = sin[ | ■ (x + (— 1)* +1 y) ], respectively. We mention that for the 
case of real coupling constants in the Hamiltonian and p = p\ = p<± ^ q = q\ = qi 
presented here, the system of equations Eqs.(|9|,|iCi|) admits solutions only for — stf tX+y > 0, 
sign[tf tX +y] = sign[tf t y- x ] and t 2 x = t 2 y . From Eq.(|63"D it can be seen that there exist at 
least one k value, for which = (for example, k x = k y = for s < 0, or k y = 0, k x = n 
for s > ). So indeed, the norm of the ground-state wave function is vanishing, it is not 
possible to remove this property as in p = q case, and as a consequence, the solution is not 
proper for describing the presented situation. 

As a conclusion, we can see that for pure real coupling constants in the Hamiltonian, the 
described solution at 3/4 filling represents a completely localized state. The part of the wave 



^ j i iuh \^\\\::"- iJuuizA-iw jjuu ijli ihiumlilvi;^ wii cvv j \ 
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function | ) that gives nonzero norm introduces on every site of the lattice rigorously the 



same number of electrons. Given by this reason, a hopping q ff Cj j(T , fl a /j j(J , or c\ a /j )(T , with 
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i 7^ j creates a state orthogonal to | ^ ) . As a consequence, all ground-state expectation 
values connected to the movement of particles within the system (i.e. kinetic energy terms 
and neighboring hybridizations) are zero. In this respect, the obtained state represents a 
paramagnetic Mott insulator, and the ground-state energy can be obtained as a sum of the 
ground-state expectation values of the on-site terms from the Hamiltonian. 

Before continuing, we would like to underline the extreme sensitivity of the solutions 
contained in Eqs.(|9|JI^) to lattice distortions. In specially we stress, that in the case of the 
undistorted unit cell (i.e. t h)X = t b , y , t b>x+y = t btV - x , V hx = V hy , V 2 ,x+ y = V 2 ,y- X ) the 
system of equations Eqs.(|| [14]) admits only pure real solutions for p and q, i.e. a completely 
localized ground-state. 



VI. THE NON-LOCALIZED SOLUTION 

In order to obtain another type of solution than that presented in the previous Section, 
we must consider the hybridization coupling constants imaginary, the Hamiltonian remaining 
hermitian. We further consider in this Section all hopping matrix elements real and VJ:* = 
V^ c for all hybridization matrix elements 7. Again, as in the case of the completely localized 
solution, two cases emerge, namely \p\ 7^ \q\, and \p\ = \q\, which will be analyzed separately. 



A. The possible solutions for imaginary p and q 

We study now the solutions allowed by Eqs. fl9j|l4D in case of imaginary p and q. The first 
group of solutions obtained, emerge at \p\ 7^ \q\, case that will be denoted as non-symmetric 
below. 



1. The non- symmetric case 



Introducing the notations from Eq.fjS^) and being interested only in situations described 
by Eq. (|54D , for |p| 7^ \q\ the solutions allowed by Eqs.(|9|,|l4l) are characterized by five indepen- 
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dent and free starting parameters, namely V2, x + y = iV2, x + y , ^2, y -x = iV2,y- x , tf, x + y , tf, y - x , and 
U, where V 1 are pure real variables. Solutions are obtained for \tf jX \ = \tf >y \ = 2^tf )X+y tf >y - x , 
9 = -tf, x + y tf )X /tf >y > and gives V ljX = (p+q*)t f)X /2, V hy = (p+q)t f>y /2, V = 0. We obtain 
a l f = xe**, a 2 j = 2x(t fty -. x /tf ty )e i<t ', a 3J = x(tf, x /tf,y)e^, a 4>/ = 2x(t fty - x /t ftX )e^, where 
X = V0, and <\> is an arbitrary phase. From this, K = 2 \x\ 2 ( \p\ 2 + \q\ 2 \2tf ty ^ x /tf ty \ 2 ), Kj = 
K — Ef = 2 \x\ 2 ( 1 + \2 tf ty _ x /tf :y \ 2 ). Introducing the notations r = tf tX /( 2 tf tX+y ), 
s = tf,y/tf, x -, for |X c (/c)| 2 we obtain 

\X c (k)\ 2 = 4 \ P \ 2 \t f , x+y \ (c^g + S ~f^C 2> _ s ^ , (64) 

where C i t has been introduced in Eq . (|63D . Since for a well defined non-zero norm we have 
to have |X c (/c)| 2 > for all k (see for example Eq. (|69D , or the explication presented after 
Eq.(f49"D), Eq.(^) shows that the non-symmetric case described here fails to represent a 
proper physical solution. 

2. The symmetric case 

From mathematical point of view the symmetric \p\ = \q\ solution is more complicated 
than the non-symmetric one. For this case two situations emerge, namely p = q, and p = q*, 
respectively. These two situations are however physically equivalent, and can be obtained 
each from other by a rotation with tt/2 of the system of coordinates. Because of this reason, 
we have to analyze in detail only one of them, namely the p = q, p = —p*. As for the 
non-symmetric solutions, we have Vy = i Vy for all hybridization matrix elements. 

The p = q solution emerge only for distorted unit cell. Five parameters can be inde- 
pendently chosen, namely tf sX+y , tf jV - x , V2, x + y , Vq, and U, so that tf }X+y ^ tf,y-x and 
sign(tf >x+y ) = sign(tf )y _ x ). The solution gives via the p = V2 >x + y /tf )X+y parameter 
the relations t CyX+y = p 2 t f)X+y , t c , y - x = p 2 tf, y - x , t c>x = \p\ 2 t f>x , t c>y = -\p\ 2 t f>y , 
Vi :X = 0, Vi tV = p tf jV , V% v -x = P tf, y -x, t 2 x = 4 tfe+y tf >y - x , and tf :V = 
-tf, x [sign(t ftX+y )} [1 + (9) 2 } 1/2 so that - //.. r //.,,// /.,-. y > 0, -r t f>y > 0, where 
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r = tf, x /(2tf jX+y ), = V tf jX+ y/[2V2, x+ y {tf jV - x - tf }X+y )]. Forthea;,/ coefficients we 
find |aij| 2 — w — u, |a 2 ,/| 2 = — r 2 (w — u), \a^j\ 2 — — w — u, |a 4 j| 2 = r 2 (w — u), 
where w = Vq/[2 p ( r 2 — 1 ) ] and u = £/ j2/ /( 2 r ). We have to consider \p\ 2 > 1 and 
r 2 7^ 1, the solution being present in the {U, Ef, Vo} parameter space on the surface^ 

U + E f = ^(M 2 - l)(r 2 + 1). (65) 



Introducing the notation s = — tf jX+y /\tf }X+y \, for z = tf yX /tf )V we obtain z = s/yl + (#) 2 , 
so < |^| < 1 emerge. For |X C (A;)| 2 , using the notation T>^ = (C 1 _ s g + r s C 2 _ s ^ ) 2 , we 
find 



|X C (A;)| 2 = 4|p| 2 |^ +y | 



1 - 




z\ 


2 


z 





v k ^ rr~ [ 1 + r<2 ~ 2 r cos (K 



(66) 



For \r\ ^ 1 the expression from Eq.(|66|) is strictly positive, so the norm of the ground-state 
wave function is nonzero and well defined. In this case we obtain besides |X c (/c)| 2 > as 
well > 0, since |X c (fc)| 2 and \Xf(k)\ 2 change their sign at the same k value. As explained 
in Sec. III., in the present case A > 0, which physically means that the diagonalized bands 
from Eq. (P0|) are non- intersecting and completely separated. 

From mathematical point of view, at the level of |X c (/c)| 2 , the main difference between 
the non-symmetric and symmetric case presented in this Section arises from the fact that 
for p ^ q we have |a l c | = |a 3)C |, |a 2)C | = |o4, c |j an d the phase of ai jC , a 3jC , (or a 2iC , a4, c ), 



differs only by or tc. Such properties are not present for p — q. In view of Eq . (|3^) , this 
means that in the non-symmetric case, coupling together (a' l c , a' 3 c ), and (a' 2 , a' 4 ) in the 
expression of |X C (A;)| 2 = K — e|, the trigonometric factors C i s g automatically arise, leading 
to Eq.(^). For the symmetric case, a such type of direct term grouping in |X c (/c)| 2 is no 
more possible. 



B. Ground-state expectation values 



Once in conditions presented in Sec. VI. A. 2. the norm of the ground-state wave function 
is well defined and |X c (fc)| 2 > 0, Jig Ag > 0, we have find a completely different solution 
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in comparison with the ground-state presented in Sec. V. Being motivated by the interest 
to find the physical properties of the system in the analyzed case, we can start now the 
calculation of all ground-state expectation values of interest. Since the study can now be 
easier done in the k coordinates, we are presenting the calculations using the k representation 
of the ground-state wave function and Hamiltonian described in Sec. III. B. 

The deduction of the main expectation values based on which the physical interpretations 
are made can be followed using the presented Appendix in their full generality in the case 
of arbitrary a^ a coefficients. However, in order to be easier for the reader to follow the 
main steps of the deduction, we are presenting below a simplified version of the calculation 
obtained in the case of site independent a- ha coefficients, which leads to the same results. 
In fact, the Appendix shows for example, that rigorously, in order to calculate ground-state 
expectation values of different Hamiltonian terms (i.e. expectation values summed over the 
spin index), is enough to consider site independence for the coefficients entering in the 
operator present in the ground-state wave function. We mention, that from physical 
point of view, the difference between a ijCT = a a , and site dependent ai i(J , is that the first 
case describes the maximal total spin S part of the ground-state wave function only, and 
the second case takes into account the full expression of the ground-state wave function. 

Introducing for the a-^ = a a case the notation \a\ 2 = I a ° \ 2 , the operator 
becomes F^ = Oi where 




(67) 



the coefficients a a being arbitrary. In the k space the contribution of F^ (taking into 
account normalized factors) becomes in this case 

F(3) = pi 3) , w here 



t? = ES fl • (68) 



and for the norm of the ground-state wave function from Eq.(fj9"|) we obtain 




(69) 



k 
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In the case of 3/4 filling, the desired expectation values becomes (we have here (...) 



\a\' 


l \X c {k)\ 2 + 


a a \* 


l \Xf{k)\ 2 




a 


2 ^ k 



Id 2 A 





! |*/(£)| 2 + 




! «£)| 2 




a: 





</L%.) = -S#- (70) 



A' 

Effectuating the sum over the spin index in Eq. (|70|) , as expected from the Appendix, the 
numerical coefficients a a disappear from the expectation values 

t 2|X c (fc)| 2 + |X / (fc)| 2 ~ , N _ 2|X / (fc)| 2 +|X c (fc)| 2 

^fc cr A; 



1} 



£<4%, CT > = -^- ( 71 ) 



Comparing with Eqs.( [A12] , |A13[ ), we see that the expectation values from Eq.([7T|) are true for 



arbitrary aii i0 ., so are correct also in the case of using the general Eq. (|20|) instead of Eq.(|T 
into p( 3 \ The ground-state expectation values of different Hamiltonian terms become 

it) = £ JtP l A '=(*)l 2 + \X;$)\ 2 ) ■ (Tf) = £ (6 ^~ £/) (2 \x f (k)\ 2 + \x c {k)\ 2 ) . 

jfe * fc * 

To have more insight about the physical behavior of the system, all ground-state expectation 

values at U > relevant for our study can be explicitly expressed from Eq. ( [72]) . The analysis 
of the described state first of all shows that N = J2k a [ ( ( fl fa a + % o- ) ) 1 = ^ as 
must be for 3/4 filling from Eq. (|7T|) ■ Secondly, summing up all contributions from Eq. ([72]) , 
we obtain for the ground-state energy E g = ( ) [ e £ ( 2 |^c(^)| 2 + |^/(^)| 2 ) + 
el(2\X f (k)\ 2 + |X c (fc)| 2 ) - 2|V^| 2 - U(\X c (k)\ 2 + \X f (k)\ 2 )]. Taking into account 
from Eq.(H) A^ = ( |X C (£;)| 2 + |X/(A;)| 2 ), this expression becomes E g = 2 Eg ( e| + 
e£ ) - Eg ( A^ 1 ) [ e| IX^)! 2 + ej. \X c (k)\* - 2 |V^| 2 ] - C7 N A . The first sum in E g 
gives 2 £y A/a (we note at this step that at U ^ we have Ef instead of Ef in the relation 
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e{ = Et + el following Eq.(E^)). Using Eqs.(|l,ffn,0), the second sum becomes — K 
As a consequence, the ground-state energy obtained from Eq . (|72|) is E g /N^ = 2Ef + U — K, 
which is exactly the value Eq /N\ given by Eq.flTSD for = 3 N\. 

Deducing now from Eq.(|71|) the total momentum distribution function = J2b,a (^| )> 
where b = c, / , we obtain = 3, i.e. a completely uniform and continuous behavior 
for the whole A;-space. This value of contains however also the contribution of the lower 
band, which being completely filled up, is 2 for every (k, a) state. This can be seen as 
well from Eq.( [All| ) of the Appendix, taking into consideration that ft f% a + ct c% a = 

kia B k,i,T + B k2a B k,2,^ writm § tne general ground state as | ) = J2er ^cr | ), 
where acr are arbitrary coefficients, taking into consideration that as shown in Eq . (]45|) , Bl 
creates a particle in the lower band, and obtaining ( J2<r B^ 2cr ) = 2 independent 

on the (k, a) indices. As a consequence, = 3 means in fact, that for the upper band 
(denoted by 1) we have n~ = 1 independent on k. This information is also contained in 
Eq.( [A10| ) of the Appendix, which gives ( Bl B^ la ) = 1, where as shown through 
Eq.(^6|), Bl creates a particle in the upper band. Using Eqs. ([All| , |A15| ) we have even 
(B% 1(T Bfixa) = l/2> ( B%2a k2a ) = 1; which via Eq.(pET|) can be transformed into the 
n| and nl momentum distribution functions defined for the starting operators / and c. 



We note at this step, that even the simplified calculation presented in Eq. flT0|) gives back 
all essential features presented above since it leads to (Bl 2 ^B^ 2a ) = 1, ( J2a Bj: t B% 1 a ) = 
1, ( Bl la B %l a ) = const. The same results are re-obtained in case of k dependent at 
coefficients in Eq.fl68l). 

Since > in the studied case, and all individual contributions in listed above 
are continuous together with their derivatives of any order in the whole momentum space. 
As a consequence, the system is in a non-Fermi liquid (NFL) state. From physical point of 
view, this property is clearly given by the presence at U > of a partially filled completely 
flat upper band, which is not hybridized with the lower band that contains dispersion. 

In order to further analyze the described state three integrals must be introduced 
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\X f (k)\ 2 \X f (k)\ 2 \X c (k)\ 2 v % 
/o - ^ ' Jl - ^ Ar ' h - \ At ■ (73) 



Using Eq.([7|), from Eq.(|72D we now find 

( T c > = -KI + h, (T> ) = — {K — Ef — U ) N A + (K — Ef — U) I + Ix , 
(V) = -2I X + (U */ 2 + VqZJ ), (£) = f/J , 

( ) = - ( V * h + V P 2 ), (E f ) = E f N A + E f I , (74) 

where, starting from Eq.([73|) Jo > 0, I\ > 0, and based on Eqs.(|| 0) we have K > 0, 
( K — Ef — U ) > 0, respectively. Introducing as in Eq . (|61~D the contribution of the 
on-site, and hopping-type Hamiltonian terms into the ground-state energy, we get ( R[ oc ) = 
E f N A + Efl - (V *I 2 + V I$), and (R mov ) = (E f -K)N A - E f I + (V * I 2 + V I 2 ) 
from where, as expected, Eq = ( Ri oc ) + ( R mov ) = (2 Ef + U — K) N A is re-obtained. 
From here one can write 

( Rloc ) = Eq + J7" , ( Rmov ) = yJ i 

J = ( K - E f - U ) N A + ( E f + U ) I - (V * I 2 + V I 2 ) . (75) 

The study of J shows that J > holds. In order to see this, via A^ = |X c (/c)| 2 + \Xf(k)\ 2 , 
we obtain J = £g Pg/Ag, where P z = K (K - el) + (K-E f )(K- e|) - {V£V n + c.c). 
Using now Eq.©, we get P % = {K - E f )\X c (k)\ 2 + K\X f (k)\ 2 + [X c (k) X}(k) V * + c.c}. 
Based on this relation, introducing the notations di iC = a i)C Xj(k), d+j = caj X*(k), and 
using Eqs.(P,|14|) we find 

4 4 

P k = H [\ d iJ? + \ d iA 2 ~ d i,cd*j - dijd* c = \ d i,c - dij | 2 > 0. (76) 

i=l ' i=l 

From Eq.(^), one can see that for all k, the value of P^ is non-negative, i.e. J > 00. So 
we have 

( R loc ) > Eg, ( R mov ) < . (77) 

Taking into consideration Eq.([T7]), it can be seen that the ground-state energy cannot be 
expressed as a sum over expectation values of on-site contribution terms of the Hamiltonian. 
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We have ( Ri oc ) > , and as a consequence, the system is not localized. From the other 
side, since (R mov ) < 0, the sum of the expectation values of Hamiltonian terms that preserve 
the movement of the particles within the system, is nonzero and negative. As a consequence, 
the ground-state energy is exactly , because in this way the system is maintaining its 
itinerant character that allows to reach the state with the minimum possible energy. In 
this conditions, the ground-state is a 2D normal state (i.e. non-symmetry broken) NFL, 
which is paramagnetic and non-insulating. The presented ground-state expectation values 
are correct only for U > 0, since for the noninteracting U = case, the | ) contained in 
Eq.flrSD represents only a negligible fraction from the linear combination of wave functions 
that build up the ground-state | ) defined at the end of Sec. II. via the arbitrary operator 
Q instead of F^ 3 \ As a consequence, we can clearly state, that the interacting ground-state 
cannot be obtained perturbatively from the ground-state of the U — case. This confirms 
the general belief, that a NFL emergence in normal phase and two dimensions has to be a 
completely non-perturbative effect, similar to ID caseil. 

Doping the system above 3/4 filling, the ground-state wave function becomes | *&g d) from 
Eq.(H). With the operator F^ in \^ d ), given by the product Il^=i present in the 
ground-state wave function, we further have [G + UP'] | ) = 0. The band remaining 
flat and being not possible to arrange the particles in such a way, to have the same number 
of electrons on every site of the lattice, the system is not localized, and remains a non-Fermi 
liquid as well. 



C. Excited states 

Starting from Eq.([EJ), it can be seen that excited states are obtained by removing -B| 2 ct 
operators from the first product of Eq . (|49|) , changing their band index to 1 (i.e. removing 
particles from the completely filled lower band to the partially filled upper band), and 
leaving intact the F^ component of the wave function. For example, removing one B-t 2 , 
we obtain 
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n %„ 

k 



II 



( ^ ) ^ (3) I ) . (78) 



Since in Eq. (f78"D the operator remains intact, it further preserves at least one /-electron 
on every site of the lattice, so P' | ) = also in this case. Using Eq. (|48|) , the remaining 
part of the Hamiltonian besides P' is essentially Y,% & ( l % i h l <j + t-x k2<y)i 

for which, the component of | ) from Eq.(|78"D which is orthogonal to the ground-state 

rci 

represents an eigenstate. Other eigenstates of this type can be obtained removing more 
than one B? operator from the lower band and changing their band index to 1. A similar 

k,2,a 

procedure can be applied also in the presence of doping. These excited states (for example, 
for one removed B>t ) , have energy E% greater than the ground-state energy with a value 
of order A = Min[ Ar]. Since in the analyzed case A^ > 0, the obtained energy spectrum 
for Eq.([T£|) is gaped. 

The excited states of the type presented in Eq. (|78|) were obtained by modifications in 
the first two operatorial components F^ and of the ground-state wave function from 
Eq . (|18D . In principle, is not possible to exclude excited states obtained by modifications 
made at the level of F^ part of the wave function, by introducing in it a double occupancy 
- empty site pair in direct space (creating in fact a supplementary double occupied site at 
the level of | ^ ) from Eq.(|l^).). In this case, the decompositions used for H presented 
in Eqs. (0,0) are not usable, and even the band structure described by E^ { , i = 1,2 in 
Eq.(pOl) is questionable at U ^ 0. However, these excited states have energy greater than 
the ground-state energy with a value of order 0(U), since a supplementary double occupied 
state has been introduced in the direct space on a given lattice site. Taking 

U » A , (79) 

the low lying excitation spectrum of the system will be clearly dominated by given by 

ki 

the states presented in Eq.(|78|). This means, that the low lying excitation spectrum will be 
gaped, so will be clearly visible in the physical properties of the system at T > 0. The 
gap becomes A^ (i.e. is k dependent), its minimum value being A. The gap symmetry 
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is a possible symmetry allowed by the described 2D lattice and depends on the starting 
parameters of the system via Eqs. (|25| . |4T| ) . 

VII. SUMMARY AND CONCLUSIONS 

In conditions in which even the exact solution for the ID periodic Anderson model 
(PAM) is not known at finite U, we present in this paper rigorous and exact solutions for 
2D - PAM in the interacting case. The described solutions are present on two surfaces of the 
T = parameter space of the model. Both solutions describe the interacting U > model, 
and the deduced ground-state wave functions cannot be obtained perturbatively from the 
non-interacting case. 

The first solution described in Sec.V. emerges for pure real hybridization coupling con- 
stants at 3/4 filling and represents a paramagnetic Mott insulator, the ground-state being 
completely localized. 

The second solution presented in Sec. VI. represents a new non-Fermi liquid in 2D 
normal (non-symmetry broken) phase and emerges in case of pure imaginary hybridization 
coupling constants (the Hamiltonian remaining hermitian) at N/N\ > 3/4 filling. This 
phase presents an n% momentum distribution function which is continuous together with its 
derivatives of any order, has a well defined Fermi energy, but the Fermi momentum is not 
definable and the system has no Fermi surface in the fc-space. At high on-site repulsion U, 
the described phase has a gap in the density of low lying excitations, and in the parameter 
space it emerges in the vicinity of a Mott insulating phase. The ground-state presents a 
large spin degeneracy and is paramagnetic. The behavior physically is given by a partially 
filled upper completely flat band situated above a normal band with dispersion. Low lying 
excitations do not give quasi-particles above the Fermi level, and have the effect to increase 
the number of particles at the Fermi energy. In this process particles are removed from the 
lower band and introduced into the upper flat band. 

We have to underline that in the case of non-distorted unit cell, from the presented 
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ground-state solutions only the completely localized state emerges. Introducing distortions 
in the unit cell (maintaining however its translation invariance) we allow in fact the emer- 
gence of the described non-Fermi liquid state in the normal phase. In this way, the presented 
non- Fermi liquid phase is strongly related to the presence of lattice distortions. 
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FIGURES CAPTIONS 

Fig.l. The hopping amplitudes for c-electrons (A), and hybridization coupling constants 
for c^f type hybridization terms (B) connecting the nearest and next-nearest neighbors 
lattice sites. The hopping amplitudes for f-electrons, and hybridization couplings for yc 
type hybridization terms are similar. 

Fig. 2. The effect of the Hubbard interaction U in the flattening of the E^ la band. The 
U and Ej^ 1 values are in |£/, x +j/| units, and k x , k y cover the first Brillouin zone. For the 
values of the other Hamiltonian parameters see text. 
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APPENDIX A: 



This Appendix contains the mathematical details related to the calculation of ground- 
state expectation values in case of arbitrary, and site dependent coefficients a- ha introduced 
by the operator P^ in Eq.(|T7]). The calculations are presented for the non-doped case. 

First of all, we compute the one-particle ground state expectation values of the form 
©6,6' = Ha b~ where b , b' = c , /, and we show that these are independent from the 

coefficients aj jt7 . 

For shake of simplicity of the presentation we use the notation cr = (cr ; )^\, <j\ £ {Til}) 
and denote by the F^ operator containing a concrete set of values = 6 au i and 
a i,l = &ai,i (i-e. zero or one depending on cr). We prove that the set of vectors defined by 
all possible cr values (see also Eq.(^9|)) 



^ ) = pW F& p® | ) 



n(%, T %a 



/J0|o> 



(Al) 



give the same expectation values for 0^/ independent on cr. We mention that the wave 



vectors presented in Eq . (|Al|) span the subspace of the ground state, they are non-orthogonal, 
but as seen below are linearly independent. 

We express also the operator F^ in terms of Bl . . For this reason first the Fourier 
transforms of the original bi 7<7 operators, then the Eq.([|7]) must be used. Due to the first 
term in Eq.( |Al|) every B*- 2 from P^ cancels out and we find 



} 



n (%, T b{ 2A 



n(%, T %i 



n f E e 



A 



0) 



Yl s s*(k),s*(cr) Z( k ) Y( k, a ) F ki i 



) . (A2) 



Here k denotes a set of k values whose first iV-j- elements (i.e. fcj for spin-up states) and also 



the last A 



1 



A A 



Aj elements (i.e. ki for spin-down states) are ordered, but there is no 



any relation between spin-up and spin-down k values. As a consequence k = { ( ki < k% < 
. . . < k N} ) , ( A;at t+ i < k N ^ +2 < . . . < k NA ) }, and £) k means a sum over all different { k } 
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configurations allowed by cr (i.e. S z ( k ) = ( — JVj, ) / 2 and S' z ( cr ) = J2i \ ( ^<ti,t — ^cr ; ,j, ) 
are equal). Denoting iVf = ATa / 2 + S z ( k ) we have 

^ = n n b\ z ( k) = n^M, ( a 3) 

- -- - ,=14/ iii 

2 t<j 2 



iV A 

f ( k, cr ) = ( - 1 £ ( - 1 ) |P| exp( i £ fc w f Q( , } ) . (A4) 

PeG i=i 

In Eq. ( |A4j ) |Q| is the parity of the permutation Q producing the spin-up spin-down separa- 
tion (i.e. Q{a) = ( er Q(i ) )£ A 1 = { T> T> •••> t> !, I, •••>!} )> an d \P\ is the parity of the 
permutation P G G = S^i,...,^} x 5{Ar T +i,...,jv A } which rearranges the k subsets j and j. 
On the other hand, the Fourier transform of i^ 3 -* = fli //„. gives 

^ 3) i o ) = e ww*) y ( k ^ ) n 4 >T n 4 4 1 o > , (as) 



which excepting Z( k ) is exactly the second parentheses from Eq. (|A2|) . Based on this 
observation it can be proved that the non-orthogonal set of ground-state wave-vectors from 
Eq. flAT]) are linearly independent. In order to do this, we prove that there exists a dual basis 
(i.e., | CT ) are linearly independent). The dual basis has the form 



k 



E <W),s*«r) n k, <r ) F k ,i | ) , (A6) 



and the proof via Eq . (|A5|) is simple 
/ vir^ 0- ' I \& u ) = 

\ 9 I 9,tr I 

= E E <W),SV) F *( k '> ) <W),*W ^( k ) F( k, CT ) ( | 4 ;1 F kll I ) 

k' k z l K J 

= E ^ z (k),S-(o-') $S«(k),fl*(aO ^*( k, Cr' ) F( k, CT ) = E &S»(o),S'{a') Y *( k , <r' ) V*( k, cr 
k k 

= E E ^(k'),5-(<r') y*( k', cr' ) 5s*(K),s*{<r) Y(k,cr) x 

k' k 



<°i( n 4,t n 4J f ( n 4 T n /iji°> 

= (0\F^F^\0) = 5^ (A7) 
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Due to the fact that all F kl | ) are eigenvectors of the operator Bl^^B^ X(j with the same 
eigenvalue than the ]l j <i^ +5z(k) f\ ^ n i> ^ +g z (k) f\ A I ) eigenvectors of the operator 
ft f^^ similarly to Eq . (|A7|) we find 

( ^ I E B,^ I ^ ) = ( I Fj« £ ^ h,„ ^ | ) 

(7 (7 

= ^ E e-^^ E ( I ^ f & /j> ^ 3) I o ) 

iYA jj' 

= E e-^-V <Wj *W <W = <W • (A8) 
Na jj' 

The dual-basis to basis expectation value from Eq. (|A5|) gives however 



/ $ u I V fit \-$r u ) = 

a 

E ( I K<r> > ( *f I E ^\ i CT V I ) = ( <° I *SU ) > ( A9 ) 



from where we find that independent from a we have 

/ \b u I V Hi \^ u ) 



1 . (A10) 



\ g,tr I ^9,<t / 

After this step, using the form of the ground state wave function from Eq. ( |A1|) we simply 
obtain for all a and cr 

Bl 5,-9 = l*!L), ^1 |*D=0. (All) 



From Eq.qSig), ( ^ | | ) / ( \ ^ > 1 , and 



%i<r I = relations arise independent on cr spin arrangement of 

the ground-state. Eq.( |All|) shows that in order to calculate from Q b b i the desired expec- 
tation values ( | Ea S'- | ) / ( #g >a . | ), first of all we may use Eq. 



to transform the initial operators into Bt (j = 1,2). After this step the non-diagonal 
Bt l B^ 2a terms cancel out based on Eq. (|All|) , and the remaining parts of the expectation 
value can be deduced from Eq.( |A10| ) and first relation of Eq.( [All| ). For example in the case 
of ct cr we find 

k,a k,cr 
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fc \ g,cr I g,cr I a 

+ X f (k)X)(k) B\ Xc B Sfha + X:(k)X}(k) Bl 2a B %lu + 

+x c (^)x / (^)4 1(T 4>,.) 



= 2|x c (fc)| 2 + ix^fc)! 2 = 2|x c (fc)| 2 + ix^fc)! 2 = x g - 6| 

|X c (fc)| 2 + iX^fc)! 2 
where we used Eq.(P^) and Eq.(^T|). Similarly 

(V u \ T ft fr I V u ) K - f f 

— 1 + 



(A12) 



(A13) 



With Eqs. ( |A12|JA13|) the independence of the analyzed expectation values on a i)(7 has been 



demonstrated. 

We concentrate now on ground-state expectation values from operatorial terms not 
summed over the spin index a. In conditions in which | ) build up a non-orthogonalized 
basis, the best way for this is to consider the T — > limit of temperature dependent expec- 
tation values. We have 

Km Tr A = Tl ^ A = ( g I A | ) 
T ™o Tre-pd Tr Wo l E CT ( ^ i ^ ) 

= ^E(^ ,CT |^I^)- (A14) 
Based on Eq. (|A14j) the ground-state expectation value of Si £>g 1 ct becomes 



(Si B- ) = — — V ( ^ Lr '°' I fit 5- I ^ ) 

\ £,1,(7 k,l,a / 2 N A ^ K 9 \ k,l,a k,l,a \ * g,v / 

= ^E(oi^ )t 44^i 3) lo) 



2^ iV A - 

_L £ e-^-V) tfjj, ^_ £ = \ (A15) 
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where the validity of the second equality is preserved by the same argument as used for the 
first equality from Eq. (|A8|) , and the last equality holds because half of the all <x has crj = a 
value. 

We mention, that the presented expectation values for Bt B^- operators are valid 
also for the complete ground-state wave function | ^ ) = J^a- 7o I )> where 7^ are 
numerical coefficients. For the case of Eq.( |All|) this is trivial, and in the case of Eq. (|A14|) 
automatically all contributions | are taken into account. For the the case of Eq.( [A10| ), 

using Eq . (|A8|) , we have 

/f^l V B~t B- I y u ) = 

a 

E jUi^ ( *,V 1 ) ( l E^L. %v I ) = 
E 7:17^(^1 1^3)^2 = 

E 7^7^ ( ^ I *^ > = K I ) , (A16) 

o\o 2 

from where ( ^ I £ CT Bt Br, I ^ ) / ( ^ I ^ ) = 1 arises. 
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